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ABSTRACT 

With a handful of measurements of limb-darkening coefficients, galactic microlensing has 
already proven to be a powerful technique for studying atmospheres of distant stars. Survey 
campaigns such as OGLE-III are capable of providing ~ 10 suitable target stars per year 
that undergo microlensing events involving passages over the caustic created by a binary lens, 
which last from a few hours to a few days and allow to resolve the stellar atmosphere by 
frequent broadband photometry. For a caustic exit lasting 12 h and a photometric precision 
of 1 .5 %, a moderate sampling interval of 30 min (corresponding to ~ 25-30 data points) 
is sufficient for providing a reliable measurement of the linear limb-darkening coefficient T 
with an uncertainty of ~ 8%, which reduces to ~ 3 % for a reduced sampling interval of 
6 min for the surroundings of the end of the caustic exit. While some additional points over 
the remaining parts of the lightcurve are highly valuable, a denser sampling in these regions 
provides little improvement. Unless an accuracy of less than 5 % is desired, limb-darkening 
coefficients for several filters can be obtained or observing time can be spent on other targets 
during the same night. The adoption of an inappropriate stellar brightness profile as well 
as the effect of acceleration between source and caustic yield distinguishable characteristic 
systematics in the model residuals. Acceleration effects are unlikely to affect the lightcurve 
significantly for most events, although a free acceleration parameter blurs the limb-darkening 
measurement if the passage duration cannot be accurately determined. 
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1 INTRODUCTION 

The light we receive from a star originates from layers at different 
depths from its surface, where the contribution of inner layers de- 
creases from the center towards the limb yielding the phenomenon 
of limb darkening of the observed surface brightness across the stel- 
lar face, which reflects the temperature of the stellar atmosphere as 
function of distance from the center. Decreasing temperatures from 
the center to the surface cause a stronger effect for shorter wave- 
lengths, making stars bluer near their center and redder near their 
limb. 

Using caustics of gravitational lenses to measure the bright- 
ness profile of closely-aligned background sources was proposed 
bv lSchneider & W eifi 1 1987) with the view toward active galactic 
nuclei, and several authors have later discussed its application to 
source stars undergoing galactic microlensi ng events that involve 
fold-caustic passages (e.g. iRhie & Bennet3ll999l : iGaudi & Gouldl 
1 999 ). Making use of both the large total magnification and a strong 
differential magnification across the stellar face during the caus- 
tic passage, galactic microlensing has successfully been demon- 
strated to be a powerful technique for revealing the brightness pro- 
file of distant stars by providing several measurements of limb- 
darkening coefficients characterizing linear or square-root laws 



lAlbrow et alJll999all20odl200lUAfonso et alJl200(tlFields et alJ 
l2003tlYoo et all2003l) with uncertaintes of a few percent. 

As discussed in detail by Dominik] 120031) . the lightcurve of 
a source star in the vicinity of a fold caustic is completely char- 
acterized by local properties, so that its brightness profile can be 
obtained already from data taken around the caustic passage with- 
out the need to find a model for the full lightcurve, which describes 
all properties of the lens system and involves a larger number of pa- 
rameters along with possible parameter ambiguities. Nevertheless, 
data outside the caustic passage region provide stronger constraints 
on the differential blending between observing sites and filters as 
well as a better assessment of acceleration effects between source 
and caustic, which can be caused by the parallactic motion of the 
Earth around the Sun or by orbital motion of the binary lens or a 
possible binary source. This additional information can be used to 
reduce uncertainties in the measurement of the stellar brightness 
profile and its characterizing limb-darkening coefficients. 

Some aspects of the measurement of limb-darkening coef- 
ficients from microlensing events involving fold-caustic passages 
have previously been studied bv lRhie & Bennettl ll999t) . This arti- 
cle provides a deeper and furthergoing discussion on this issue and 
on some general properties of lightcurves of microlensed sources 
near fold caustics with respect to revealing their brightness profiles. 
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Sect.[2]addresses the relation between the lightcurve and the 
stellar brightness profile and discusses the sensitivity of microlens- 
ing as function of phase during the caustic passage and fractional 
radius of the source star. With linear limb darkening as example, 
Sect. [3] deals with parametrized profiles. A simulation of typical 
sampled data for different sampling rates is described in Sect. [4] 
which are used in Sect.|5|to investigate the amount of information 
about the linear limb-darkening coefficient contained in different 
regions of the lightcurve. In Sect. [6] consequences arising from 
the adoption of an inadequate limb-darkening profile profile are 
discussed, whereas Sect. Q addresses the influence of an effective 
acceleration between source and caustic. The paper finishes with 
a summary of the results and recommendations on the observing 
strategy in Sect.|8| 




(3) 



(4) 



and K denotes the complete elliptical integral of first kind (c.f. 

Consider now a star with angular radius p*#E whose leading 
limb enters or whose trailing limb exits a fold caustic at time t* t and 
which needs the timespan 2 1+ to cross. The parameter 



v = ±- 



t-tf 



+ 1 



(5) 



2 INTENSITY PROFILE AND LIGHTCURVE 

Let us consider a gravitational lens of total mass M at distance Dl 
from the observer acting on a source at distance Ds- A characteris- 
tic scale is then given by the angular Einstein radius 



4GM D s - -Dl 



D L D s 



(1) 



while its corresponding physical radius at the lens distance is the 
Einstein radius te = -Dl 9e- 

With the combined magnification due to the two critical im- 
ages of a point-like source at a distance u±9e from a fold caustic 
given by 



^ rit («x)=(^) V2 ©(«±) 



(2) 



where Rt describes the caustic strength, 111 > (u± < 0) for a 
source inside (outside) the caustic, and O(x) is the step function, 
the corresponding magnification for a circle of angular radius v8e 
with its center at u±8e reads 



then describes the phase of the caustic passage, where the upper 
sign corresponds to a caustic entry and the lower sign corresponds 
to a caustic exit. The source is completely outside the caustic for 
77 < 0, completely inside the caustic for 77 > 2, and passes the 
caustic for 77 2, where rj — 1 corresponds to the center of the 
star crossing the caustic. As will be shown in the next section, the 
reference time t{ is preferred over the time tf = tj ± when the 
source crosses the caustic, because it can be immediately identified 
from a distinctive feature in the lightcurve. 

With p denoting the fractional stellar radius, the stellar bright- 
ness profile for a given filter reads 



(6) 



where / denotes the average surface brightness and is a di- 
mensionless function describing the radial profile, so that 



(7) 



The magnification of the critical images takes the form 
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A crit (r?;£)= ( ^ ) Gf (»?;£), 



1/2 



(8) 



and the characteristic caustic profile function Gf (77; £) can be writ- 
ten as 



(9) 



where 1"(ry, p) describes the response produced by the caustic to 
the brightness profile £(p). In general, the determination of the 
brightness profile £(p) involves an inversion of Eq. which con- 
stitutes a Fredholm integral equation of first kind. The solution of 
this problem is non-trivial and has been discussed in some detail by 
iHevrovskd 120031 . 
Since 



A cri t(^;e) = 2 / A3t(V—L,pp* ) Z(p)pdp, (10) 



one finds that 



T(r),p)=2p 



1/2 , 



n-i 



(ii) 



Eq. 1101 reflects the fact that, in order to account for the fraction 
of light provided by each circle, the surface brightness at the corre- 
sponding fractional radius has to be weighted by its circumference, 
which is proportional to p. 

The function T(rj, p) is shown in Fig.Qas a function of the 
caustic passage phase 77 for selected values of the fractional stel- 
lar radius p, and as a function of the fractional stellar radius p for 
selected values of the caustic passage phase 77. As long as the cir- 
cle corresponding to the fractional radius p is outside the caustic, 
only its non-critical images contribute to its light, so that T(r], p) 
is zero for for 77 < rj° ut — 1 — p, where the caustic is touched and 
the sensitivity jumps to the value (2p) 1//2 . The sensitivity increases 
with 77 for ri° ut < 77 < rff until the circle is completely inside the 
caustic at rfp = 1 + p, where T(rj, p) diverges. For 77 > rfp , the 
sensitivity asymptotically approaches zero as T(r], p) ~ pr)^ 1 ^ 2 . 

The sensitivity of the lightcurve to the brightness profile at 
different fractional radii p for a given caustic passage phase 77 qual- 
itatively depends on whether the souce is completely outside the 
caustic (77 < 0), completely inside the caustic (77 = 2), or in the 
process of passing it with its center outside (0 ^ 77 < 1), inside 
(1 < f] ^ 2), or at the caustic (77 = 1). For any caustic pas- 
sage phase 77, the sensitivity T(r\, p) vanishes at the source cen- 
ter (p = 0), while it vanishes identically for any fractional radius 
p if the source is completely outside the caustic (77 < 0), where 
none of the corresponding circles have critical images. With parts 
of the source moving inside while the source center is still out- 
side the caustic (77 < 1), only the outer parts of the source with 
p > p^ ut — 1 — rj have non-vanishing sensitivities. T(rj, p) jumps 
from zero to [2(1 — 77)] 1 / 2 at the fractional radius p^ ut and increases 
towards larger p. For the source center being at the caustic (77 = 1), 
one finds that T(l, p) oc p 1/2 . If the source moves further inside 
(1 < 77 ^ 2), increasing fractional radii subsequently dominate 
the sensitivity function through their infinite peak at pJJ 1 = r) — 1, 
to which T(t), p) rises monotonically for p < pjj 1 . For p > p™, 
T(r), p) shows a monotonic decrease from the infinite peak larger 
77, while a local minimum exists for smaller rj, so that a decrease is 
followed by an increase towards p — 1. Finally, for sources com- 
pletely inside the caustic (77 > 2), T(rj, p) rises monotonically 



with the fractional radius p from zero at the source center (p = 0) 
to a maximum at the limb (p — 1) and becomes asymptotically 
proportional to p for large 77, for which T(rj, p) — prj~ 1/l2 . 

The properties of the brightness profile sensitivity function 
T(j), p) show that photometric observations during the passage of 
a source star over a fold caustic provide a one-dimensional scan 
of its brightness profile £(p), during which each fractional radius 
p is most efficiently probed in different ways on the two occa- 
sions where its associated circle touches the caustic, namely at 
Vp Ut — 1 — P' where the profile sensitivity function jumps from 
zero to the finite value (2p) 1//2 , and at = 1 + p, where it be- 
comes infinite. Therefore, outer parts of the source are most effi- 
ciently probed near the beginning or the end of the caustic passage, 
whereas inner parts are most efficiently probed while the center 
of the source passes the caustic. Influenced by the fact that circles 
corresponding to larger fractional radii have a larger circumference 
and therefore contribute more light for the same surface brightness, 
the integrated sensitivity of the caustic to the brightness profile over 
the full course of the caustic passage increases with fractional ra- 
dius, so that data of constant quality would yield more information 
about the outer than about the inner parts of the source. 1 

Rather than Gf (77; £) itself, one observes the lightcurve (e.g. 
lDominikl2603l) 



TTlfold(t) = 7T7f 



2.5 lg < 1 + 



fit 



+ Gather (±Wf (t 




a + 



(12) 



in the vicinity of the caustic passage, with t£ and tj" as defined 
above, t being an arbitrarily chosen unit time, and raf the observed 
magnitude at time t* { . With ms ^ m\ denoting the intrinsic source 
magnitude, the parameter g\ reads 



1/2 



1/2 



(13) 



Finally, Gather describes the temporal variation of the magnifica- 
tion due to images of the source star that do not become critical if it 
hits the caustic, which is measured by o> f * and can be approximated 

as 



Gothcr(A,<7) = g (exp{A/(?} - 1) 
~ A (|A/ff|«l). 



(14) 



The establishment of Gf (rj; £), from which the stellar bright- 
ness profile £(p) is to be inferred, therefore requires a determina- 



tion of the 5 parameters tf , t t 



(; f , and i2> t from a fit to the 



collected data, where in particular, t f and t+ determine the relation 
between the caustic passage phase 77 and the elapsed time t. 



3 PARAMETRIZED PROFILES 

As an example of a parametrized stellar brightness profile, let us 
consider the case of linear limb darkening characterized by the co- 
efficient ^ r < 1, where 



For constant exposure time during the caustic passage, larger magnifica- 
tions however account for a higher photometric precision. 
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Figure 2. Stellar brightness profiles £(p;T) = (1— IT) § {0 } (p)4T (p) 
for linear limb darkening as a function of the fractional stellar radius p, 
where the two extreme cases of uniform brightness and maximal limb dark- 
ening, corresponding to T = or V = 1, respectively, are shown. 



£(p;r) = (i-rK {0} (p) + r£ {1} ( P ), 

with the basis functions 



(15) 

(16) 
(17) 



While £{o}(p) corresponds to a uniformly bright source, £{i}(p) 
is linear in cos i) = y/l — p 2 , where i? denotes the emergent an- 
gle. These basis brightness profile functions are displayed in Fig.|2| 
Any profile £(p; F) of the form given by Eq. I15i is a linear super- 
position of these curves and falls between them. 

For such brightness profiles, the caustic profile function reads 



GHn; r) = (l - r) G}, {0} ( v ) + rG?, {1} (i/) , 



(18) 



where Gf { y{v) an d @t {lyiw) denote the caustic profile func- 
tions for the basis stellar brightness profile functions, which can 
be expressed by means of analytic or semi-analytic functions as 
( Schn eider & Weilj Il987t ISchneider & WaeoneJ Il987l: iDominikl 



V < 

Hp.' 

for < 77 < 2 
^[(2-^(^27^) - 




(19) 



-(l-r,)E [j2/r, 



for r) 5? 2 



where K and E denote the complete elliptical integrals of the first 
or second kind, and 



G 



f,{i} 



fa) 



f (5 - 277) »r 

§ [(5- 
+ (1 



3/2 

f 277) (V 



+ 



2) 



3/2] 



for 77 < 
for < 77 < 2 

for 77 > 2 



(20) 



respectively, and are shown in Fig.|3|together with their derivatives 
Gf ' = dGf /dr] and Gf " = d 2 Gf /dr/ 2 . For any choice of T, 
Gf (77; r) vanishes for 77 ^ 0, and rises to a peak, which is shifted 
towards smaller 77 for stronger limb darkening. Asymptotically, the 
decrease after the peak behaves as Gf(?7; r) ~ (77— 1)~ 1,/2 regard- 
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Figure 3. The caustic profile function Gf (77; T) = (1 — T) G£ ^ y(v) + 
r G* r x i (77) and its derivatives Gf '(77; T) and Gf "(77; T) as a function of 
the caustic passage phase 77. Solid lines correspond to uniform brightness 
(r = 0) and dashed lines to maximal limb darkening (r = 1), while 
general profiles with < V < 1 fall between these two extreme case. 



less of the adopted stellar brightness profile, which corresponds to 
the behaviour of a point-like source. While the slope of Gf tniv) 
vanishes for 77 = and remains finite for 77 = 2, it shows a jump 
discontinuity at 77 = and a negative infinite value at 77 = 2 for all 
caustic profiles Gf (77; F) with r / 1. The trend of initially smaller 
values of the profile function Gf (17; V) and its slope for stronger 
limb darkening is reversed on the rise to the peak. The slope shows 
strong variations for small 77 > and near 77 = 2, and little vari- 
ation in the remaining regions. For 77 = 0, the change of slope 
with 77 becomes infinite unless F = 0, for which a finite negative 
value is taken. Approaching 77 = 2 from larger values, the slope of 
Gf (77; r) diverges to positive infinity, whereas it diverges to nega- 
tive infinity for the approach from smaller values, except for F = 1, 
where a finite negative value is approached. 

A change in the linear limb-darkening coefficient dF causes 
the brightness profile £(p; F) to vary at each fractional radius 
9 by [<9£(p;F)/dr]dr. Due to the linearity of £(77; V) in F, 
9£(p; F)/dF equals the difference between the basis profiles, i.e. 
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Figure 4. Difference between profile functions AGf (77) = Gf — 
Gj ^qj(^) and its derivatives AG^'(rj) and AGj"(r;) measuring the sen- 
sitivity of the lightcurve to the limb-darkening coefficient T as a function of 
the caustic passage phase r?. 



ggfoT) 

ar 



e { i}(p)-?{o } (p) = Ac(p) : 



(21) 



while all higher derivatives of £(p; T) vanish identically. 2 As Fig. [2] 
shows, A£(p) decreases monotonically from A£(0) = 1/2 to 
A£(l) = -1. With dA£(p)/dp = -(3/2) p (1 - p 2 )~ 3/2 , its 
slope becomes flat at the source center and decreases monotoni- 
cally towards the limb where it becomes infinite. The normalization 
forces A£(p) to change it sign, which occurs for po = v5/3 w 
0.745. 

Changes in £(p; V) caused by variations of V alter the caustic 
profile function Gf (77; F) according to Eq. ^9} with the sensitivity 
of the caustic profile at a given passage phase 77 to these varia- 
tions characterized by T{rj, p) as discussed in the previous sec- 
tion. Gf (77; r) inherits the linearity in V from the brightness profile 
C(p;r),sothat 



dGUv,T) 
dT 



Gf, { i } (r/)-G f * {0 }(»7) = AG? fa). 



(22) 



measures the sensitivity of the lightcurve to variations of the limb- 
darkening coefficient V and for its derivatives, one finds in analogy 
dGt'(mr)/dT = AG?' (77) and dG$"(m Y)/dV = AG?"(rj). 
These functions are displayed in Fig.|4] 

Since the stellar magnitude mfdd(t) is a function of elapsed 
time t rather than of caustic passage phase 77, a reliable measure- 
ment of the limb-darkening coefficient V requires not only a suf- 
ficient sensitivity of the caustic profile function Gf (77; T) to the 
variation of T, but in addition also the determination of the param- 
eters if , , TOf , fff , and o>f , which establish the relation between 
Gf (77 : r) and mf id(/j), as pointed out in the previous section. 



2 Such a property does not hold for the popular parametrization of linear 
limb darkening where £(p; c) = (1 — c/3) — 1 ^1 — c + c 1 — p 2 ^ . 



The determination of these parameters from the photometric data 
over a given region of the lightcurve can be viewed in relation to 
the local assessment of the caustic profile function Gf (77; T) and 
its derivatives Gf '(77; F) and Gf "(77; Y) with respect to the caus- 
tic passage phase 77 as well as their derivatives with respect to T, 
namely AGf (77), AGf '(77), and AGf "(77), where the amount of 
information provided increases with the absolute value of each of 
these quantities. 

With T(r),p) ^ 0, the sign change in A£ at po = VE/3 
causes AGf (77) to change its sign at 77 ( ° ut = 0.716 and ?7 n = 
1.762, where AGf (77) < for < 77 < 77g ut or 77 > 77 n , whereas 
AGf (77) > for 77§ ut < 77 < 770 1 . A change in V causes the largest 
relative changes in the observed flux at the extrema of AGf , which 
occur at <? n = 0.263, 77 max = 1.276, and 77^ = 1.971, i.e. 
close to the point at which the caustic is touched by the limb of the 
source located inside. 

By means of the sensitivity function T(rj,p), variations of 
77 during the caustic passage alter the weight given to the differ- 
ent fractional radii p and AGf '(77)] increases with the absolute 
variation in A£(p) for the specific radii that dominate T(rj, p) for 
a given caustic passage phase 77. With dA£(p)/dpj Af (p) < 
decreasing monotonically from the souce center to the limb, and 
moreover, the caustic producing a weaker response characterized 
by T(rj, p) to variations of the brightness profile near the source 
center, the variations in AGf (77) are small for 77 ~ 1 and large 
as 77 — > or 77 — > 2. In fact, Fig. [4] shows that | AGf '(77)! as- 
sumes large values for small positive 77 and in the close vicinity of 
77 ~ 2, and becomes infinite at 77 = and 77 = 2, whereas much 
smaller values are assumed for the remaining part of the caustic 
passage, and the variation of AGf '(77), given by AGf "(77), be- 
haves in a similar way. Roughly speaking, when comparing differ- 
ent phases of the caustic passage, the largest response to a variation 
of the limb-darkening coefficient is caused around the time when 
the stellar limb touches the caustic with the source center being 
inside (77 ~ 2) with all AGf (77), AGf '(77), and AGf "(77) assum- 
ing reasonably large absolute values or even tending to infinity as 
77 — > 2. Next comes the region over which the stellar limb touches 
the caustic while the source center is outside, where in contrast 
AGf '(77)1 decreases to zero as positive 77 — > and vanishes iden- 
tically for 77 $; 0. In contrast, much smaller effects result while the 
inner parts of the source pass the caustic, for which both AGf '(77) 
and AGf "(77) are close to zero or j AGf (77)! is small if one of the 
former quantities tends to assume a reasonable non-zero value. 

In order to assess the potential of different regions of an 
observed lightcurve for the measurement of the (linear) limb- 
darkening coefficient T, one needs to consider not only the flux 
variation with F at the corresponding caustic phase, but also the 
observed flux and its uncertainty. Let F{ denote the flux at time t{ , 
so that 7nfoi d (t) = ra\ — — 2.5 ]g[Ffoid(t) / Fif]. The ratio ^(77) 
between signal-to-noise for flux deviations caused by a variation 
of the limb-darkening coefficient V and signal-to-noise for the flux 
Ff is given by 



z( v ) = 



dF fc 



(t?±r,t±;T) 



or 



F*a 



(23) 



With the assumption that the flux uncertainties a follow Poisson 
statistics, so that a = a\ [Fj Ff) 1 ^ 2 , where a\ is the uncertainty 
of Ff, and neglecting the temporal flux variation in non-critical 
images (i.e. adopting u)f = 0), one obtains 



Z{rj) = ar 



AGf (77) 



y/a T G*(t7;F) + 1 



(24) 
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Figure 5. Ratio Z{rf) between signal-to-noise for the flux variation 
[c9.F(ij + rj t^- ; T) / 8T] / a caused by altering the limb-darkening co- 
efficient r and signal-to-noise Ff/ai for the flux measurement out- 
side the caustic for selected caustic-peak-to-outside magnification ratios 



peak 



where 



IF*. 



or = 



(Fpeak/Ff*) 



G 



(25) 



pcak,r 



and Gpcak.r is the maximum of G} (77; V) which corresponds to the 
observed flux F pca k- 

Using r = 1 as reference, Z(rj) is shown in Fig.[5]for selected 
values of the peak-to-outside flux ratio Fp^/Ff. Compared to the 
effects on the flux difference as displayed in Fig. [4] caustic pas- 
sage phases near the beginning of an entry or the end of an exit are 
favoured, yielding a maximal signal for phases close to - but not at 
- this point. However, [AGf '(rj)\ is small in this region, whereas 
it peaks for r\ — 2, i.e. when the caustic entry is completed or the 
caustic exit starts. 

For Z(rj) = 20 (a* /Ff), a single data point implies an abso- 
lute uncertainty of 0.05 for the limb-darkening coefficient T. If one 
considers F pcs k/Ff ~ 8 as for recently observed events, Fig. [5] 
yields \Z(rf)\ ~ 0.25, so that for a photometric precision of a few 
percent at the caustic outside, tens of data points yield a reliable 
measurement of Y. 

The identification and assessment of some distinctive charac- 
teristic local features in the light curve that do not depend on the 
stellar brightness profile allow the precise determination of indi- 
vidual lightcurve parameters, leading to a reduction in parameter 
degeneracies and uncertainties. For any stellar brightness profile 
that does not vanish at the stellar limb, the slope of the lightcurve 
exhibits a jump discontinuity at r\ — and its identification in the 
lightcurve directly yields m* and t* . Another distinctive feature for 
such brightness profiles is the infinite slope and sign change of the 
curvature at r\ — 2, where temporal identification of both r\ — 
and rj = 2 yield ft- With Gf (??;r) = for 77 < 0, data taken 
while the source is completely outside the caustic provide a)f, and 
with m* and g* can then be measured from the asymptotic form 
G*(rj; F) ~ (77 — l) -1 / 2 for sources inside the caustic. The possi- 
bility of measuring each of the parameters t*, m*, and tlif directly 
from data obtained around the passage of the stellar limb over the 
caustic while the source is outside (77 ~ 0) makes this phase of the 
caustic passage more valuable than that of the limb passage while 
the source is inside (77 ~ 2), which alone does not provide any di- 
rect measurements of individual light curve parameters. Only the 
combination of both of these phases yields direct measurements of 
tt and g* . In contrast, data taken while the inner parts of the source 
pass the caustic (77 ~ 1) do not provide any direct measurements of 




Figure 6. Simulated dataset for the lightcurve parameters t* = 0, tj- = 
6h,m* = 16.3, g* = 0.05 (i = 1 h), u>* = -0.001 h -1 , and T = 0.5. 
The standard sampling interval is At = 6 min, while data corresponding 
to a sparser sampling with At = 30 min only are displayed in bold. The 
measurement uncertainty has been assumed to be 0.015 mag for and 
scaled to other magnitudes assuming the flux to follow Poisson statistics. 
More details can be found in the text. The division of the dataset into an 
outside region A, a central region B, with subregions Bi and B2, and an 
inside region C, is indicated by vertical dashed lines. 



lightcurve parameters neither alone nor in combination with data in 
other regions. 



4 SIMULATED PHOTOMETRIC DATA 

4.1 Creation of dataset 

In order to study the uncertainties of the measurement of the lin- 
ear limb-darkening coefficient in dependence of the sampling of 
the lightcurve for different regions, a simulated dataset for a typi- 
cal caustic exit in accordance with the capabilities of the PLANET 
campaign 1 Domi nik et alj200 2) has been created. Contrary to caus- 
tic entries, which have to be caught by chance, caustic exits can 
be predicted and therefore provide regular opportunities to obtain 
high-quality data during caustic passages. For a G/K subgiant in the 
Galactic bulge at a distance Ds ~ 8.5 kpc with a brightness I ~ 17 
and a radius R ~ 3 7?*, the angular radius becomes 8* ~ 2 /ias, 
so that the choice of t* = 6 h corresponds to a proper motion 
perpendicular to the caustic of (j, ~ 7 /iasd -1 and to a velocity 
v ~ 70 km s _1 at the lens distance Dl ~ 6 kpc. For the unit time 
t — 1 h and the adopted limb-darkening coefficient F — 0.5, the 
choice = 0.05 yields a relative magnification ^ pea k/^4f ~ 12.5 
of the caustic peak relative to the caustic exit (at time t* ), for which 
we assume rrif = 16.3. Finally, u>* — —0.001 h _1 and t* = are 
adopted. 

The sampling is characterized by the sampling interval At = 
6 min, an uncertainty a At = 45 s of the point of time at which 
the measurement is taken, and the sampling phase shift t p hasc = 

1.2 min. The uncertainties in the observed magnitudes have been 
assumed to be 00 = 0.015 at ra\ (i.e. the relative uncertainty of the 
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measured flux is ~ 1.5 %) and scaled to other magnitudes, so that 
the uncertainty for the corresponding flux follows Poisson statistics 
in being proportional to the square-root of its value. Moreover, the 
uncertainty of the magnitude at given time has been smeared with 
a relative standard deviation of f„ — 0.125. 

For tmi n — —15 h and t max = imin + nAt = 1.2 h, a 



synthetic dataset (t 



a^') containing n — 162 simulated 



measurements of the time of observation, the magnitude, and its 
uncertainty has thus been created using 



,(0 



U{t 

fa ao 



<T() 



! + i At + tphasc, CAt) • 
10 {0.2[m*-m fold (i<*))]} 



10{°' 2 [ m f- m fold(* <l) )]} 



m w = AA(m fold (t W ),a W ) 



(26) 

(27) 
(28) 



where Af(fi, cr) denotes a value drawn randomly from a normal 
distribution with mean [i and standard deviation a. 

In order to simulate possible losses in the data acquisition due 
to observing conditions or technical problems, data points have 
been removed from the dataset at random with a probability of 
qi = 5%. By random selection of data points with a probability 
of q s = 20 %, a restricted dataset has been obtained, which simu- 
lates a sparser sampling with an expected average sampling interval 
of 30 min. The data used for the subsequent analysis are shown in 
Fig. |6| with 27 points, which belong to the sparse sampling only 
shown in bold, while the dense sampling includes all 153 points 
shown in bold and in light. 

Missing or sparse coverage over part of a caustic passage are 
likely to occur if its duration exceeds the time of coverage possi- 
ble with a single telescope. Although the coverage by more than 
one telescope would require different parameters m* f ' s ' and 
(c.f. Dominik 2003), a careful assessment of parallax between the 
sites <AnetalJl2 002). and different achievable photometric preci- 
sions, this will be neglected in the following for reasons of simplic- 
ity, where site-independent fit parameters m\ and g\ will be used, 
which is equivalent to treating the data as resulting from observa- 
tions with a single telescope. 



4.2 Tests on statistical properties 

In order to provide a reference and to check the created data set, let 
us have a look some of its statistics with respect to the underlying 
set of true parame ters p tl .„„. 

Let us follow lRhie & Bennett! Il999h in using the caustic pas- 
sage phases at which AG\(r\) vanishes, namely ry{J ut = 0.716 and 
7/ n — 1.762 as boundaries to divide the data into an outside region 
A (V < r lo nt )> a central region B (r/g ut < 77 < r^ 11 ), and an inside 
region C (77 Vo 1 )- Furthermore, let us subdivide the central region 
B at 77S = (77g ut + r ?0 ,1 )/2 into two regions Bi (^ ut < 77 < r?g) 
and B2 (r)o ^ rj < t?™)- These regions are indicated in Fig. [6] 

TableQgives the number of data points N for the full dataset 
and each of the regions A, Bi, B2, and C for both samplings as 
well as Xtvuc = X 2 (Ptruo). the ratio Xttue/N, which is expected 
to be unity, and the probability P(x 2 ^ Xtrue)- While the standard 
upper-case letters refer to the simulated data for the dense sam- 
pling with At — 6 min, the sparse sampling for a certain region is 
indicated by the corresponding lower-case letter. One sees that the 
average absolute size of the residuals differs among the different 
regions, in particular region B2 shows rather large residuals while 
the adjoining region C partly compensates this with smaller than 



Table 1. Number of data points TV and scatter of residuals for different 
regions of the simulated dataset. 



Region 


-^truc 


N 


Xl uc /N 




A 


51.5 


50 


1.02 


0.42 


Bi 


34.5 


31 


1.11 


0.30 


B_> 


42.8 


30 


1.43 


0.06 


C 


36.7 


42 


0.87 


0.67 


ABC 


165.5 


153 


1.08 


0.23 


a 


5.69 


11 


0.52 


0.89 


hi 


3.40 


3 


1.13 


0.33 




4.38 


5 


0.88 


0.50 


c 


2.70 


8 


0.34 


0.95 


abc 


16.17 


27 


0.60 


0.95 



X 2 (ptrue) ls me vame of x 2 for me true parameters p true and 

= p(x 2 > 



Xtrue)- Uppper-case letters refer to the dense and 



lower-case letters to the sparse sampling. 



expected residuals. At first sight, it might be surprising to find a 
probability of only about 6% for region Bi. However, this means 
that the probability that at least one of four regions shows a prob- 
ability of that order or less is ~ 25%, so that about a quarter of 
all simulated datasets should share this feature, which is reflected 
in the probability for the full dataset ABC. The fraction of points 
selected for the sparse sampling is ~ 18 %, somewhat smaller than 
the expected 20 %. In particular, only ~ 10 % of the points in re- 
gion B2 made it to b2. Regions a and c show residuals whose ab- 
solute values are on average much smaller than expected, with a 
probability of being smaller than this between 5 and 10 %. Two re- 
gions with probabilities of about 10 % make the sparse sample less 
likely than the dense one with P(x 2 ^ Xtrue) ~ 5 %. 

Systematic deviations in the data can be quantitatively as- 
sessed by means of a run test, which is based on the statistics of the 
signs of the residuals irrespective of their absolute values, whereas 
a x 2 test is based on the absolute values while being blind to the 
signs. Let a 'run' being defined as the longest contiguous sequence 
of residuals with the same sign, and let N denote the total number 
of data points, iV+ the number of points with positive residuals, and 
iV_ the number of points with negative residuals. For N > 10, the 
distribution of the number of runs n r can be fairly approximated by 
a normal distribution with the expectation value 

(29) 



£{n T ) = 1 + 



N 



and the standard deviation 



2N+N- 



N 



N 2 (N- 



(30) 



Results of run tests on the residuals of data in all regions for both 
adopted samplings with respect to the true lightcurve are listed in 
Table|2| which do not contradict the assumption that the scatter in 
the residuals is random. 

Without knowledge of the true model parameters, the assess- 
ment of the goodness-of-fit by means of \ 2 is hampered by the 
fact that the reported error bars frequently do not properly reflect 
the true error bars of the measurement (e.g. lUdalski et alJI 19941 
lAlbrow et aT]|l998h forcing the introduction of modelled scaling 
factors and/or syste matic errors fe. g I Ansarill 1 994t |Dominikll 1 996t 
Tsaoras et al. 2003), so that the hypothesis of underestimated error 
bars rather stands against rejecting the model on a basis of a larger 
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Table 2. Run tests on the lightcurve residuals for the true model parameters. 



Region 


N 


N+ 


N— 


£(n r ) 


cr(n r ) 


n° hs 


S 




ABC 


153 


68 


85 


76.56 


6.09 


75 


0.26 


0.40 


abc 


27 


14 


13 


14.48 


2.54 


12 


0.98 


0.16 



From the total number of N data points, N+ show positive and TV— show negative residuals, implying £(n T ) runs with a standard deviation <r(n r ) to be 
expected. With n* s observed runs, the deviation S = [£(n T ) — n° hB ]/a(n T ) is related to a probability P r = P(n r $C n° ). 



than expected Xmin- This is not a problem for the assessment of the 
run test, which does not depend on the size of the error bars. 



5 MEASUREMENT OF PARAMETERS 

Without going into much detail. lRhie & Bennetj fl999^ have stated 
that a reliable measurement of the limb-darkening coefficient re- 
quires a reasonable coverage of the lightcurve in each of the regions 
A, B, and C. Using the dataset presented in the previous subsec- 
tions, let us study the information content of each of these regions 
for the two different realizable sampling rates and the adopted pho- 
tometric precision that can be expected from the PLANET cam- 
paign. 

Table|3|gives the parameters obtained by the minimization of 
X 2 for the simulated data shown in Fig.|S|in selected regions. The 
quoted parameter uncertainties 8p refer to projections of the hy- 
persurface A\ 2 = x 2 (Pmin + S P) - X 2 (P m in) = 1 onto the Pa- 
rameter axes, where p min denotes the parameters at the x 2 min- 
imum Xmin — X 2 (Pmin)' which correspond to intervals contain- 
ing a probability of 68.3 %. Since the model is not linear in the 
parameters, upper and lower bounds differ, and 95.4 %-intervals 
(Ax 2 = 4) are not twice the size of 68.3 %-intervals. Table |4] 
shows the corresponding values of Xmin> me probability P x 2 — 
P(x 2 ^ Xmin) indicating the goodness-of-fit, and the x 2 excess 
Axt 2 ruc = X 2 (Ptmc) ^ X 2 (Pmin) of the true parameters p truo with 
respect to the best-fit parameters p min . The values of Axtrue have 
to be seen with regard to Xtrue as listed in Table Q for the differ- 
ent regions of the lightcurve and the average misestimation of error 
bars indicated by Xtme/N. 

Densely sampled data (every 6 min) for all of the regions 
(ABC) provide a precise measurement of the limb-darkening coef- 
ficient r with an uncertainty of ~ 2.5 % and for a sparser sampling 
(every 30 min) over all regions (abc), the uncertainty still remains 
below ~ 6 %. For the dense (sparse) sampling, the uncertainties on 
the time tf marking the end of the caustic exit and on the passage 
half-duration t^r are ~ 30 s (~ 90 s), while the flux at tf is deter- 
mined within ~ 0.4 % (~ 1.5 %) and the relative uncertainties on 
gi and d) f * turn out to be ~ 0.5 % (~ 2 %) and ~ 10 % (~ 40 %), 
respectively. 

From data in the outside region A alone, T is determined with 
an uncertainty of ~ 8 %, While the slope discontinuity at the end 
of the caustic exit and the data thereafter provide reliable measure- 
ments of tf, m\, and (b\ with uncertainties not that much in ex- 
cess of those for the full dataset, the uncertainties in g\ and 
are largely increased by factors of ~ 20 or ~ 60, respectively. By 
combining the outside region A with the inside region C, the latter 
parameters are accurately determined and the uncertainties on all 
parameters are close to those for the full dataset (ABC), the uncer- 
tainty in r being ~ 4 %, so that the central region B does not add 
much valuable information if data in both of the regions A and C 
are available. 



Data taken in each of the inside region C or the central re- 
gion B alone is insufficient to provide a reliable measurement of 
the limb-darkening coefficient, where the central region B makes 
the worse case. However, already a combination of the inside re- 
gion C with the adjoining half of the central region B2 yields V 
with an uncertainty of ~ 10 %, which improves to ~ 6 % if the full 
region B is included. The gain in the parameter uncertainties of the 
other parameters by adding region B2 to ABi is rather moderate, 
and the uncertainty on the prediction of the caustic exit time tf from 
data in regions ABi is ~ 30 min, while the uncertainty on the cor- 
responding magnitude is ~ 2 mag. Compared to the case in which 
only data in region A are used, the uncertainties on tf, m*, g f * and 
uif are at least an order of magnitude larger, while the uncertainty 
on t+ however is about 3 times smaller. 

Although data in region AC provide equal or better constraints 
than data in region AB or ABi for all of the parameters tf , t* , mf, 
<7f and u)f , region ABi provides a measurement on V of similar 
quality as region AC, and region AB even provides the same small 
uncertainty of ~ 2.5 % as the full dataset ABC, making data in 
region C obsolete in this case. Adding region Bi to A allows a reli- 
able determination of t+ , and significantly reduces the uncertainty 
in gj? , whereas the uncertainty in tf shows little and those in rrif 
and u>f show no significant improvement. The main gain from re- 
gion B 2 is in reducing the uncertainty in tj", whereas rrif and (if 
as well as their uncertainties remain practically unchanged, and the 
improvements on tf and g\* are quite small to negligible. 

With dense sampling of region A yielding an uncertainty of 
~ 8 % on F and regions B or C being unable to provide a reason- 
able measurement, none of the regions alone could yield a sufficient 
accuracy with the sparse sampling. From the combinations of two 
of these regions, ab yields a marginal limb-darkening measurement 
with an uncertainty of ~ 20 %, while the uncertainty is ~ 15 % 
for be and ~ 10 % for ac. Nevertheless, one sees that it is possible 
to obtain a measurement of the linear limb-darkening coefficient V 
already with a small number of data points with high photometric 
accuracy over the course of the caustic passage. 

The comparison of fits with at most one of the regions A, B, 
or C being densely and the remaining parts of the source sparsely 
sampled shows that only dense sampling of region A yields a rea- 
sonable reduction of the uncertainty of all parameters, where the 
uncertainty of the limb-darkening coefficient r of ~ 3 % is close 
to that obtained for the densely sampled full dataset. In contrast, 
the improvements by increasing the sampling of regions B or C are 
marginal, although the better sampling of region C yields signifi- 
cant reductions in the uncertainties of tj", # f * and &f . 

As for the dense sampling of the outside region A, combining 
sparsely sampled data in this region with data in one of the other 
regions is highly valuable, where in contrast to the dense sampling, 
the inside region (C or c) provides more additional information (in 
particular through direct measurements of tjr and g}) than the cen- 
tral region (B or b) and coverage of both regions contributes to a 
significant reduction of the uncertainty on the limb-darkening co- 
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Table 3. Best-fit parameters and their uncertainties for the simulated data in selected regions of the lightcurve. 



Region t* [h] t£ [h] m* g* -&* [h" 1 ] F 



ABC 


n mm +0.0091 

-o.oooi_ 0088 




.0097 
.0095 


i a noRT+0.0040 
16.2967_ 0040 





n^n9i+ ' 00025 

• UJUZ -0.00025 


n nnnon+0.00011 
0.00099I ooon 


0. 


.5121°: 


012 
.012 


aBC 


-0.007_ 0016 


5.997_q 


.017 
OKI 


1 R OO7+0.011 

16.287_ 0011 





nc - n79 +o.ooo60 

LMU(Z_ 00060 


n nnnoo+0. 00012 
0.00099_ ooon 


0. 


.512« 


(121 
(12(1 


Abe 


-0.007_ 0010 


O.UlU_Q 


.013 
.013 


1K nnno+0 0052 
16.2998_ 0051 





04975+ 00048 

,ui»< J_ 00048 


n nmoo+0. 00032 
0.00128^ ooo32 


0. 


522+° 


016 
Dili 


aBc 


o.oo_ 002 


o.UU7_q 


.02.3 
()24 


i « onn+0-013 
16.290_ 0013 





n^n-?8+ 00082 

.UOUOO_ 00082 


n nm nn+0. 00041 
0.00120I ooo37 


0. 


■516« 


(124 
024 


abC 


n nn/i+0 020 
-0.004_ 00lg 


5.9o7_q 


.019 
(117 


1fi noo+0.011 

lb.^!S8_ Q11 





05060+ 000061 

,UOUDU_ 00060 


n nmnQ+0 00012 
0.00103_ 00012 


0. 


■518« 


025 
(125 


abc 


n nnn+0 024 

o.ooo_ 023 


o.UU7_q 


.025 
(125 


i « 9m +0.013 
16.291I 0013 





o t sn20+°' 00085 

.UDUZU_ 00085 


n nmon+0. 00045 
0.00129^ ooo40 


0. 


520lg 


031 
03(1 


ABi 


-0.003_ 0011 


o.U3o_q 


.044 
(142 


ic -jnno+0 0067 
16.3002_ 006g 





04941+ 00086 
,u<i»<i±_ 00082 


n nm ii\+0. 00051 
0.00130I ooo52 


0. 


499+g 


(12(1 
019 


AB 


n nn-j+0 010 
0.003_ 0010 


c r\r\ a tU 

O.UU4_q 


.015 
Dili 


ic -jnno+0 0067 
16.3002_ 006g 





04981+ ' 00072 

.U1DOL_ Q 00066 


n nm qq+0. 00052 
0.00133_ 0Q053 


0. 


.512« 


013 
013 


Ab 


n non+0 010 
°- 009 -0.012 


cr r»no~rO 

o.9ys_Q 


.021 
(121 


i r Qnm +0.0066 
16.3001± . 0067 





04984+ ' 00073 
.ut»oi_ 00067 


n nm qq+0. 00052 
0.00133^0.00053 


0. 


529+g 


018 
018 


aB 


n nfio+0 099 
°- 062 -0.060 


5.9»9_q 


.054 

.049 


16.41Io.il 


0422+ 0066 
U.U4ZZ_ .0052 


n nm +0.010 
O.OOIIqoos 


0. 


430lg 


(18(1 
092 


dD 


n n -+0.09 

u.u/_ 17 


r: Q7c;+0 

o.y / o_q 


.007 
(>4(i 


lb.4d_ 23 




0.0418lg;gg 8 


n m ,+0.010 
U-011_ ois 


°- 44 io. 


10 
09 


AC 


002+ 011 


6.001+g 


.011 
(110 


16 2970+ - 0040 

lo.z»(U_ 0040 


U.UOUlO_ 00027 


00102+ 00013 

u.uuiuz_ 00012 


0. 


.5171°; 


(12(1 
(12(1 


Ac 


006+ 012 

U.UUD_ 012 


6.011+p 


.013 
.014 


16.3000lg;™? 


04970+ 00049 
u.uty<u_ 00049 


00132+ 00034 
U.UUlOZ_ 00033 


0. 


.521« 


(121 
020 


aC 


-0 023+ 024 


5.9841° 


.020 
(117 


16.287tg-.gl 2 


0.0b0»^!_ 00063 


0009 I 5+ ' 00014 
u.uuuyo_ 00013 


0. 


482tg 


(13,9 

036 


ac 


-0 014+ - 028 


5.998^° 


.020 

.026 


16 291+ 013 
io.zai_ Q13 


05027+ ' 00088 

U - U0UZ '-0.(10(189 


n nni98+o ooo 47 
0.00128_ 00041 


0. 


493lg 


(144 
(142 


B 2 C 


19+ 034 
u ' iy ~0.61 




17 

31 






034+0-082 


00032+ 0045 

U.UUUOZ_ 0008 


0. 


494lg 


054 
048 


BC 


09+ 31 


6.05« 


16 
23 


16.6± 2 ; 3 




Q3O+0.060 
u.uoa_ 035 


0004+ 0032 
u.uuu^_ 0009 


0. 


.523+0 


(132 
028 


Be 


-0 31+ 70 


5.86t° 


35 
27 


K 7+21 




095+ 081 
U.Ub)O_ og5 


0.0038_ 0042 





562+g 


(152 
049 


bC 


n 25+ - 22 


6.i3±° : 


12 

24 


17 1+ 2 ' 1 
1.3 




024+0- 059 
U.UZ4_ 22 


-0.000l±g;^l 


0. 


.521« 


(138 
0.3(1 


be 


o i+°- 3 


6.08+° 


20 

.53 


16.5+ 2 ; 3 




04+°- 15 

u ' u ^-0.04 


0.0008i°;g»» 


0. 


■545+g 


063 
052 


A 


-0 010+ 016 
u.uiu_ 015 


6.66^ 


67 

49 


16 3004+ - 0068 
lO.OUU1_ oo69 


0425+ 0053 
U.U1ZO_ 0058 


00112+ - 00048 
U.UU11Z_ 00046 


0. 


466tg 


(144 
043 


B 


-0.3+i;= 


6.0+? 


8 
5 






1 5+036 
u - io _0.14 


000+ 050 
u.uuu_ 038 


(1- 


io- 6 ) 


+ 1()- 6 
-0.59 


C 


c r+l-9 
°'°-0.8 


o oq+0.94 

OOU -0.39 


1 *" 5Z -0.14 




71+0- 27 

u - ' i -0.35 


031 + 017 
U.Udl_ 018 





The simulated data are shown in Fig.lSTand corresponds to the true parameters t* = 0, tj" = 6 h, = 16.3, = 0.05 (i = 1 h), ui^ = —0.001 h 1 , and 
r = 0.5. The letter A denotes the outside region (r) < rjg ut ), B denotes the central region (r)g ut i^r] < rf^ 1 ), and C denotes the inside region (r) ^ Vq 1 )' w i m 
region B subdivided into Bi {rj^ ut s^r] < r] c ) andB2 (rj c < rj < r; n ), where AGj? vanishes at r;g ut = 0.716 andr/'^ 1 = 1.762, and rj c = (r;g ut +rf^)/2. 
Uppper-case letters refer to the dense sampling (At = 6 min), while corresponding lower-case letters refer to the sparse sampling (At = 30 min). The 
quoted uncertainties refer to intervals enclosing a probability of 68.3 %. 

efficient V. For both samplings of the outside region, most of the 
additional information contained in the other regions is already ex- 
tracted with the sparse sampling. In all cases, the sampling rate in 
region A strongly affects the accuracy to which all parameters can 
be determined. 

With data only in the inside and the central region being ob- 
served, a denser sampling of the inside region strongly reduces the 
parameter uncertainties while a denser sampling of the central re- 
gion has very little effect. 



6 INAPPROPRIATE BRIGHTNESS PROFILES 

If an inappropriate brightness profile is adopted, the residuals of the 
data with respect to the model lightcurve show characteristic sys- 
tematic trends. These are shown in Fig.Qfor models of the full sim- 
ulated dataset displayed in Fig. [5] where the linear limb-darkening 
coefficient has been fixed at F = 0.3 or F — 0.7, respectively, 
whereas r = 0.5 marks its true value. The largest residuals occur 
near the end of the caustic passage and there are long sequences 
for which all data points either show larger or smaller magnifica- 
tions than the theoretical curve. In particular, all data points near 
the time where the central parts of the source passes deviate to one 



side, whereas points where outer parts of the source pass deviate to 
the other side. The sign of the residuals is just inverted for the two 
cases r = 0.3 and F = 0.7, representing a too weak or a too strong 
linear limb-darkening term, respectively. 

The parameters and their uncertainties for these two models 
and corresponding models for the sparse sampling are shown in 
Tabled One sees that the achievement of the best-possible fit with 
the wrong brightness profile leads to an inappropriate choice of the 
caustic passage half-duration t+ and the point of time tf when the 
trailing limb of the source exits the caustic, where t+ is underes- 
timated and tf is shifted towards the caustic inside for a weaker 
limb darkening (F — 0.3), whereas t+ is overestimated and tf is 
shifted towards the caustic outside for a stronger limb darkening 
(r = 0.7). The remaining parameters m*, , and (i f * are fairly 
reproduced. 

There is no significant difference in the residuals between ei- 
ther of the two adopted samplings being applied to the different re- 
gions of data. The same systematic trends still persist even if only 
data in the outside region A are modelled, although the size of the 
residuals is decreased compared to modelling the full dataset. 

Table|6|shows the results of both \ 2 - an d run tests for the mod- 
els with fixed limb-darkening coefficient of the data corresponding 
to the two different samplings. The larger number of data points for 
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Table 5. Best-fit model parameters for the simulated dataset with fixed inappropriate limb-darkening coefficient. 



r 


Region 


t* [h] 


** M 


III * 


9* 


-£)* [h" 1 ] 


0.3 
0.3 
0.7 
0.7 


ABC 
ahc 

ABC 
abc 


-o.i38t°;°£ 

1259+ - 0046 
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Models with T = 0.3 or T = 0.7 using the full simulated dataset corresponding to dense sampling (ABC) or the reduced dataset corresponding to sparse 
sampling (abc). The quoted uncertainties correspond to 68.3 % confidence intervals. 



Table 6. \ 2 - and ran test for models with a fixed inappropriate limb-darkening coefficient. 
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Minimization of \ 2 yielded \ 2 



for the listed models, which corresponds to a probability P x 2 



data points within the specified region revealed 2V+ positive and 7V_ negative residuals forming n° bs runs, where £(n r 
deviation <r(ri r ). The discrepancy 8 = [£ (n r ) — n° bs ]/c(n r ) con'esponds to a probability P r = P(n r ^ n° bs ). 
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Figure 7. Residuals of data with respect to model lightcurves for (incorrect) 
stellar brightness profiles with V = 0.3 and T = 0.7, where the true value 
of the underlying simulated dataset (as shown in Fig.|6j is T = 0.5. 



Table 4. x 2 tes t on models based on data in different regions and corre- 
sponding x 2 excess of true parameters. 
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The probability P x 2 



P(x 2 ^ Xmi n ) provides a measure of the 



goodness-of-fit, while Ax 2 ruo = X (Ptrue) ~ X (Pmin) § ives the ex- 
cess in x 2 of the true parameter set p truc over the best-fit parameter set 

Pmin- 



the dense sampling strongly increases the significance and each of 
the tests clearly recommends rejection of both of the models with 
r = 0.3 and F = 0.7. For the sparse sampling, this the case only 
for the model with Y — 0.3, whereas the model with V — 0.7 
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This definition preserves the caustic passage taking place for ^ 
77 ^ 2 and lasting 2t* . The proper motion of the source perpen- 
dicular to the caustic reads 



t* [h] 



Figure 8. Residuals of data with respect to model lightcurves with fixed (in- 
correct) acceleration parameters aj" = 0.01 h _1 and = —0.01 h -1 , 
where the true value of the underlying simulated dataset (as shown in Fig.lol 
is &i = 0. 



might be accepted. The x 2 tests turn out to be more powerful than 
the run tests for all models. The power of the latter is limited by the 
requirement that subsequent systematic deviations to the same side 
need to be established for each of the regions over which they per- 
sist separately, while the use of the sum of the squared deviations 
in x 2 tes t s exploits the power of the total number of data points, 
which can become quite large for N > 100. However, as pointed 
out in Sect. 14.21 the uncertainty in the size of the true error bars 
limits the power of rejecting models by means of a x 2 test. 

While stellar brightness profiles described by a linear limb- 
darkening law have been studied in detail in this section, it is ob- 
vious that other profiles show a similar behaviour if the strength of 
limb darkening is under- or overestimated. 



(32) 



where (if = fJ. (tf ) and with the constant acceleration fi ± , the 
acceleration parameter a x is defined as a x — /x x //if . For a x < 
0, the source will reverse its direction of motion for sufficiently 
large t where /i ± (t) becomes zero. Requiring that the source does 
not move in opposite directions at tf ± 21* and tf implies the 
restriction a x ^ — 1/ (2t^ ) and ensures that tf ±2 t x istheearlier 
time for which rj — 2 is reached. 

For microlensing observations toward the Galactic bulge, let 
us adopt source and lens distances of Ds ~ 8.5 kpc and Z?l ~ 
6 kpc, respectively. With the Earth's velocity of ~ 30 km s _1 , a 
change of the relative proper motion perpendicular to the caustic of 
= K e D DlDs ( U e/( 1AU )) ~ 1-5 xlO -2 /iasd" 1 is induced, 
where ^ 1 is a geometrical projection factor depending on the 
angles between the Earth's acceleration vector and the line-of-sight 
and between its projection perpendicular to the line-of-sight and the 
relative proper motion between source and caustic. For a relative 
perpendicular proper motion of ji ± ~ 10 /xasd -1 , corresponding 
to a velocity of v — Dl ji ~ 100 kms~ at the position of 
the lens, the acceleration parameter becomes a x ,® = /i^ j /fi ± ~ 

4 x 10 _5 Ke h_1 - 

The orbital motion of a binary lens causes a rotation of 
its caustic with its orbital period and an alteration of shape and 
strength due to changes in the angular separation between its con- 
stituents. For a total mass M ~ 0.7 Mq and semi-major axis 
a ~ 1.6 AU, the orbital period becomes P ~ 2.4 yr, the Ein- 
stein radius for Ds ~ 8.5 kpc and Z?l ~ 6 kpc is te ~ 3.1 AU, 
and the angular Einstein radius is #e ~ 0.54 mas. For moderate 
eccentricities, the orbital acceleration does not differ much from 
that for circular orbits. Depending on the orbital inclination and 
phase, the separation between the constituents of the binary lens 
falls into the range ^ r p ^ a, where r p = a for all times 
in face-on orbits. Since a ~ 0.51 te, the lens is a close binary 
with respect to its caustic topology for any mas s ratio q, showing 
a diamond-shaped ca ustic at its center-of-mass (Erdl & Schneider 
1993; Dominik 1999). The acceleration induced by the rotation for 
a point on the caustic located at an angle 9 from the center-of-mass 
reads fi^ ot = 4rc 2 K rot (r p /a) P~ 2 9, where the factor r p /a ac- 
counts for the projection of the orbital motion onto the sky, while 
the geometrical factor « rot ^ 1 is related to the angle between the 
projected acceleration vector and the proper motion between source 
and caustic. For a mass ratio q ~ 0.8 and r D = a, the angular ex- 
tent of the caustic is ~ 0.26 #e (Dominik 1999), and for a source 
intersecting the caustic at 6 ~ 0.1 #e, the angular acceleration per- 



7 ACCELERATION EFFECTS 
7.1 Origin and order estimates 

The revolution of the Earth around the Sun and the orbital mo- 
tion of the binary lens or a possible binary source 1 Dominik 1998; 
ISmith. Mao & Paczvh ski 2003) cause an effective acceleration be- 
tween lens and caustic, which leads to a modified relation between 
the caustic passage phase rj and the elapsed time t, namely 



t-tf 



1 + a* ti 



(31) 



so that (i 1 - ~ 10 /iasd 1 implies an acceleration parameter of 
a*,rot ~ HT 5 ^ h" 1 . 

While face-on orbits provide the largest effect on the acceler- 
ation from the rotation, the alteration of the angular separation of 
the constituents is largest for edge-on orbits and causes the size 
of the caustic to change proportional to d 2 — (r p /r_B) 2 . With 
(d 2 /dt 2 )d 2 = 2(dd + d 2 ) and d(t) = (a/r E ) coa[(2n / P)t + <p ] 
for a circular edge-on orbit, the induced acceleration for a point 
on the caustic at angle 9 reads fi OBC = 8ti 2 k^ sc P~ 2 9, which be- 
comes /i t J sc = 5.4 x 10 -3 K oac /iasd _1 for the parameters cho- 
sen above, so that the acceleration parameters becomes d* osc ~ 
2 x 10~ 5 « osc h _1 , where the geometrical factor k^c ^ 1 accounts 
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Table 7. Estimates for the correlation coefficient between residuals and elapsed time. 
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For models involving the full dataset over regions A, B, and C, the estimated correlation coefficients pi and p2 correspond to data in the subregions Bi and 
B 2 , while p± = (pi ±p 2 )/2 give the antisymmetrized or symmetrized correlation for the full region B. Symmetry is indicated by a vanishing p+, while 
a vanishing p~ indicates antisymmetry. With TVi = 31 data points in the region Bi and 7V 2 = 30 data points in the region B2, the standard deviation on 
Z\ 2 = 0.5 ln[(l + pi j2 )/(l — pi 2], is cr zi £S a zi ss 0.19 and the standard deviation on = (zi ± Z2)/2 is a z ± fs: 0.13. Values that are compatible 
with zero within one standard deviation are shown in italics. 



for the orbital phase and inclination as well as for the caustic cross- 
ing angle. 

The effective acceleration caused by orbital motion strongly 
increases towards shorter binary periods, however, in order to in- 
tersect the increasingly smaller caustics, smaller impact parame- 
ters are required, which are less likely to occur. Larger acceleration 
effects are therefore not impossible but rather improbable. In gen- 
eral, the contributions of parallax and the rotational and oscillatory 
contributions by orbital motion are rougly of competitive order. 



7.2 Signature of acceleration 

In order to cause deviations of 1-2 %, the acceleration parameter 
needs to exceed |q^| > 0.01 h _1 for the parameters = 0.05 
(t — 1 h) and t* t — 6 h, which have been chosen for the simu- 
lated dataset, and is therefore required to be much larger than the 
expected values, indicating that effective acceleration is unlikely to 
play an important role. The influence of acceleration increases with 

/ acc OC (i*) 1/2 (fffT 1 

OC Q^^^ 1/2 10-°- 4(mS - m f ) 

= &i (tU^) 1/2 10-°- 4(ms - m f\ (33) 

so that stronger effects occur for longer passage half-durations t+ 
in combination with smaller angular stellar radii or smaller 
proper motions p ± , as well as for larger caustic strength Rf, al- 
though a very large increase is not likely to occur, whereas a de- 
crease results from larger ms — rrif , i.e. larger magnification at tf . 

The residuals for models with aj; = ±0.01 h -1 , as shown in 
Fig. [8] exhibit characteristic systematic trends, which can be distin- 
guished from those occuring with the adoption of an inappropriate 
stellar brightness profile. For the latter case, the residuals during 
the central part of the caustic passage (region B) are roughly sym- 
metric with respect to the boundary between the subregions Bi and 
B2 while deviations to one side dominate for adjoining parts of the 
limb regions A and C, and just the very outer limb in region A 
shows a strong deviation to the same side as the central part of re- 
gion B. In contrast, the residuals for models including an inappro- 
priate assumed acceleration are roughly antisymmetric about the 
temporal midpoint for each of the regions A, B, and C, showing a 
linear trend and reaching maxima at about the boundaries of these 
regions and the endpoints of the caustic passage. Therefore, tests of 
symmetry or antisymmetry in the different regions can be used to 



decide whether model discrepancies are due to acceleration effects 
or the adoption of inappropriate brightness profiles. 

For N data points (xi, yi), where 1 ^ i ^ N, an estimate for 
the correlation coefficient — 1 p ^ 1 is obtained as 



P = 



(34) 



(35) 



where Fisher's transform 

1 , 1 + P 
z = — m - 

2 1 - p 

asymptotically follows a normal distribution around its expectation 
value with standard deviation a z — 1/ \/ N — 3. 

For the subregions Bi and B2, estimates for the correlation 
coefficient, p\ and p 2 , and the corresponding values of Fisher's 
transform, z\ and z 2 , have been calculated for the previously dis- 
cussed models involving an incorrect brightness profile or accel- 
eration parameter based on the complete dataset ABC. The anti- 
symmetrized or symmetrized correlation over the whole region B 
is characterized by p — (pi ± p 2 )/2 and z = (zi ± Z2)/2 with 
a z ± = 0.5 y/(N! +N 2 - 6)/[(ATi - 3)(JV 2 - 3)]. For a sym- 
metry (as for an incorrect brightness profile), one expects p + to 
vanish, while p~ is expected to vanish for an antisymmetry (as 
for an incorrect acceleration parameter), whereas signatures of a 
symmetry or antisymmetry are obtained from values that signifi- 
cantly differ from zero. The obtained values are displayed in Ta- 
ble where entries in italics are compatible with zero within the 
standard deviation. The obtained estimates for the correlation coef- 
ficient reflect the expected behaviour and allow to distinguish be- 
tween the two different causes of systematic residuals unless the 
random scatter of the residuals dominates over the systematic trend, 
which is the case for |d* I = 0.001 h _1 , where both p + and p" 
are compatible with zero. 

The model parameters and their 68.3 % confidence intervals 
are shown in Table [8] One sees that larger acceleration parameters 
yield smaller limb-darkening coefficients T, shorter caustic passage 
half-durations t+ , smaller Cjf and larger g£ , whereas smaller accel- 
eration parameters alter the estimates for F, t*, Uf, and g$ in the 
opposite direction. 

While probabilities corresponding to \ 2 as shown in Table [5] 
reveal the incompatibility of the model with the data for the se- 
lected values of = ±0.01 h _1 and &i = ±0.003 h _1 due to 
the size of the absolute deviations, the results of a run test show the 
presence of systematic trends for = ±0.01 h _1 . As can be seen 
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Table 8. Best-fit model parameters obtained for selected fixed acceleration parameters. 
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Best-fit parameters and 68.3 % confidence intervals for models involving selected fixed values for the acceleration parameter ci+ using the full dataset 
corresponding to dense sampling for which a^- had been adopted. 



Table 9. x 2 - and run test for models with fixed acceleration parameters. 
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On the full set of N = 153 data points, the x 2 -minimization yielded xf n s n with P x 2 = P(x 2 X^m) and tne mn test revea l e d positive and A?_ 
negative residuals, for which £(n r ) runs with a standard deviation <r(n r ) are expected, whereas n° bs have been found, yielding 8 = [£(n r ) — n° bs ]/ff(n r ), 
corresponding to a probability P r = P{n r ^ n* s ). 



from both Table [8] and Tableland the comparison with the mod- 
els parameters and their uncertainties for vanishing acceleration as 
listed in Table [3] the effects of acceleration become negligible for 
&i = ±0.001 h" 1 . 



7.3 Influence on limb-darkening measurement 

In general, acceleration can be assessed by including as addi- 
tional free parameter in models of the lightcurve. If it cannot be 
accurately determined, degeneracies with other model parameters 
can lead to an increase of their uncertainties. Although variations of 
the amount of limb darkening and of acceleration provide different 
signatures over the full course of the caustic passage as previously 
shown, similar effects in local regions of the lightcurve can arise 
so that with a partial sampling their origin cannot be resolved and 
the measurement of the limb-darkening coefficient is blurred by the 
possible presence of significant acceleration. 

The uncertainties on best-fit model parameters including a free 
acceleration parameter ajr based on the simulated dataset for sev- 
eral selected regions of the lightcurve shown in Table fTol and the 
results of a \ 2 tes t shown in Table ITT1 demonstrate that accelera- 
tion affects models for which the caustic passage duration cannot 
be determined with sufficient accuracy, while otherwise a precise 
measurement of the acceleration leaves both the other model pa- 
rameters and their uncertainties practically unchanged and the un- 
certainty of the measurement of the linear limb-darkening coeffi- 
cient does not suffer from acceleration effects as a comparison with 
models for the same regions that neglect acceleration effects listed 
in Table|3|shows. 

Even with the adopted dense sampling of the outside region 
A alone, for which the acceleration parameter is compatible with 



values of the order of ~ 0.1 h _1 , a meaningful measurement 
of r cannot be obtained and large parameter degeneracies arise if 
acceleration is considered to be a free parameter. In contrast, if data 
in the outside region are combined with data from at least one of 
the other regions, the amount of acceleration is accurately deter- 
mined and therefore does not influence the uncertainties on other 
parameters already for the more sparse sampling. Without a dense 
sampling of the outside region, the remaining parts of the lightcurve 
need to be densely sampled in order to yield a meaningful measure- 
ment of the linear limb-darkening coefficient, where its uncertainty 
is increased by a factor of two compared to the assumption of van- 
ishing acceleration. 

Data taken over all regions allows to assess the acceleration 
parameter accurately, down to \&jr\ < 10~ 3 h _1 for the sparse 
sampling and even down to \&jr\ < 6 x 10~ 4 h _1 for the dense 
sampling, which is sufficient for acceleration effects effects on the 
lightcurve and the determination of the other parameters and their 
uncertainties to become negligible. 

In any case, parameter degeneracies involving acceleration 
can be limited by adopting an uppper limit on |A+| resulting from 
an assessment of the dynamics of the binary lens and the Earth's 
motion at the time of observations. 



8 SUMMARY AND RECOMMENDATIONS ON 
OBSERVING STRATEGY 

As laid out bv lSchneidei"& Weifi ll987l) . the differential magnifi- 
cation across the face of a source during its passage over a caus- 
tic created by an intervening gravitational lens can reveal its radial 
brightness profile through frequent accurate photometric observa- 
tions. Galactic microlensing surveys such as OGLE-III are capable 
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Table 10. Best-fit model parameters including acceleration and their uncertainties for the simulated data in selected regions. 
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The true parameters are tf = 0, t+ = 6 h, mf 
intervals enclosing a probability of 68.3 %. 



16.3, g* = 0.05 (t = 1 h), 6j* = -0.001 h _1 , a 4 - = and T = 0.5. The quoted uncertainties refer to 



Table 11. \ 2 test on models including acceleration for simulated data in 
selected regions and corresponding x 2 excess of true parameters. 
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"^•min 


d.o.f. 
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ABC 
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146 


0.20 


5.6 


abc 


11.4 


20 


0.93 
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ab 


8.4 


12 


0.75 


5.1 


ac 


5.3 


12 


0.95 


3.1 


BC 


111.8 
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0.13 


2.3 
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43.8 


43 


0.44 


7.7 



The probability P x 2 = P(x 2 X^in) P rov ides a measure of the 
goodness-of-fit, while Ax 2 ruc = X 2 (Ptruo) - X 2 (P m in) § ives the ex- 
cess in x 2 of the true parameter set p truc over the best-fit parameter set 

Pinin* 

of providing ~ 10 such passages of stars in the Galactic bulge per 
year from which a measurement of the stellar brightness profile can 
be obtained, which reflects the variation of temperature with dis- 
tance from the center providing a powerful technique for probing 
stellar atmosphere models. 

Contrary to caustic entries, the corresponding exits are pre- 
dictable from data obtained during the characteristic rise in magni- 
fication to a peak I Albrow et al. 1999b; Jaroszvnski & Mao 2001), 
which allows to schedule frequent observations before the caustic 
passage begins. 

While the caustic passages are likely to last from a few hours 
to a few days, photometric measurements with lm-class telescopes 
with an uncertainty of less than 1.5 % with sampling intervals of a 
few minutes can be obtained for stars brighter than ~ 17th mag- 
nitude as currently being carried out by the PLANET collaboration 
(Dominik et al. 2002), offering the possibility to collect a few hun- 
dred data points during the course of the passage. 

In the vicinity of the caustic passage, the lightcurve can be 
described by means of a characteristic profile function Gf (»7;£), 
which depends solely on the dimensionless normalized stellar 
brightness profile as function of the fractional radius p and the 
caustic passage phase r\. G\ (r\\ £) can be seen as the response deliv- 
ered by the caustic to the specific form of the brightness profile. The 
weight of the contribution of the brightness at a specific fractional 
radius p to the caustic profile function Gf (r/; £) is given by a func- 
tion T(r\, p), which measures the sensitivity of the lightcurve at the 
point of time that corresponds to the caustic passage phase r\ to a 
local variation of the stellar brightness profile at fractional radius 
p. An inspection of Tlj), p) shows that the caustic passage pro- 
vides a one-dimensional scan of the brightness profile where each 



fractional radius is most efficiently probed as it touches the caustic, 
which happens twice during the course of the caustic passage, so 
that the majority of information about outer radii is provided at the 
beginning and the end of the caustic passage, while inner regions of 
the source reveal their identity close to times when the source cen- 
ter passes. However, the integrated sensitivity over the full duration 
of the caustic passage increases with fractional radius. 

In general, the extraction of the stellar brightness profile £(p) 
from the observed lightcurve rri{ id(t) involves the determination 
of 5 model parameters (tf , t t , rrif, g^, ujf), which relate mf id(t) 
and the caustic profile function Gf(j)\ £), as well as the solution of 
the integral equation, which relates Gf (77; £) and £(p) by means of 
the caustic sensitivity function T{r\, p). 

The variation of coefficients characterizing a parametrized 
stellar brightness profile causes variations in the observed 
lightcurve through variations of the brightness profile at all frac- 
tional radii. The sensitivity of the lightcurve to variations in such a 
coefficient V therefore depends both on the properties of T(r), p) 
and (9£(p; r)/<9r. For linear limb darkening, this sensitivity is 
largest near the beginning of a caustic exit, followed by the sur- 
roundings of its end, whereas the smaller integrated sensitivity over 
the caustic passage for smaller fractional radii and the fact that 
d^(p;T) /dT shows little variation with p for the inner parts of 
the source imply a small sensitivity while these pass the caustic. 
However, the identification of the end of the caustic exit with the 
characteristic feature of a jump discontinuity in its slope allows a 
direct measurement of the corresponding point of time t} and mag- 
nitude m,, while data points just after the end of the caustic exit 
can be used to measure the parameter iif* directly, which charac- 
terizes the variation of magnification due to non-critical images of 
the source. For this reason, sampling of the surroundings of the end 
of the caustic exit is more valuable than of its beginning. Monitor- 
ing of both of these regions provides direct measurements of the 
caustic passage half-duration t+ and the caustic rise parameter . 

Unless a precision measurement on the linear limb-darkening 
coefficient T is desired, a sampling interval at the limit of the ca- 
pabilities of the monitoring campaign is not required. For a caustic 
passage lasting 12 h, a sampling interval of ~ 30 min (correspond- 
ing to ~ 25-30 data points) is sufficient to provide T with an uncer- 
tainty of less than ~ 8 %, while a dense sampling corresponding to 
an interval of ~6 min (with ~ 150 data points) would reduce the 
uncertainty to less than ~ 3 %. This offers the possibility for mea- 
suring limb-darkening coefficients in several broadband filters with 
the same telescope on the same microlensing event or to monitor 
other microlensing events during the same night (e.g. to look for 
anomalies caused by planets around the lens stars). 
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Since coverage of the surroundings of the end of a caustic exit 
is most effective in providing an accurate measurement of the linear 
limb-darkening coefficient, it should be the goal of any observation 
strategy to try to obtain data in this region, where the sampling 
rate should be chosen with regard to the desired accuracy of the 
measurement and the number of broadband filters. Moderate sam- 
pling over other regions of the light curve is more valuable than an 
increased sampling over the end of the caustic exit, where a very 
dense sampling over these other regions however does not provide 
much additional information but can be used to compensate for a 
missing coverage of some parts of the caustic passage. 

The adoption of an inappropriate stellar brightness profile or 
the neglect of acceleration effects, which include contributions due 
to the revolution of the Earth around the Sun or the orbital motion 
within the binary lens or a possible binary source, leads to well- 
distinguishable characteristic systematics in the model residuals. 
Although acceleration is unlikely to cause significant effects on the 
lightcurve for most events, if used as a free model parameter, how- 
ever, it blurs the measurement of the stellar brightness profile if the 
duration of the caustic passage cannot be determined with sufficient 
accuracy, which is the case if only the surroundings of the end of 
the caustic exit are sampled. Nevertheless, the related degeneracy 
can be limited by applying a reasonable upper limit to the absolute 
amount of acceleration resulting from an assessment of the dynam- 
ical properties of the binary lens and the Earth's motion at the time 
of observation. 
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